Constraints on the threshold K nuclear potential from 
FINUDA A Z(K- top ,n-)iZ spectra 

A. Cieply a , E. Friedman b , A. Gal b '*, V. Krejcink a,c 

a Nuclear Physics Institute, 25068 Rez, Czech Republic 
h Racah Institute of Physics, The Hebrew University, 91904 Jerusalem, Israel 
c Department of Physics, University of Maryland, College Park, MD 2074-2-4.111, USA 



Abstract 

Isa hypernuclear formation rates in stopped K~ reactions on several p-shell 
targets are derived from hypernuclear formation spectra measured recently 
by the FINUDA Collaboration and are compared with calculated Isa for- 
mation rates based on a chirally motivated coupled channel model. The 
calculated rates are about 15% of the derived rates, and in contrast with 
previous calculations depend weakly on the depth of the threshold K~ nu- 
clear potential. The A dependence of the calculated ls\ rates is in fair 
agreement with that of the derived ls\ rates, showing a slight preference for 
a deep density dependent potential, Re Vk-(po) ~ —(150-200) MeV, over 
a shallow potential, Re Vr- (po) ~ ~ 50 MeV. These new features originate 
from a substantial energy and density dependence found for the in-medium 
subthreshold K~n — > 7r~A branching ratio that enters the hypernuclear for- 
mation rate calculations. 
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1. Introduction 

How strong is the K~ nuclear interaction? Various scenarios proposed 
for kaon condensation in dense neutron-star matter [lj], and more recently 
for quasibound K~ nuclear clusters ji[ and for self-bound strange hadronic 
matter depend on the answer to this question which has not been resolved 



* Corresponding author: Avraham Gal, avragal@vms.huji.ac.il 



Preprint submitted to Physics Letters B 



April 8, 2011 



todate. A modern theoretical framework for the underlying low-energy KN 
interaction is provided by the leading-order Tomozawa- Weinberg vector term 
of the chiral effective Lagrangian which, in Born approximation, yields a 
moderately attractive K~ nuclear potential Vr- : 

^- = -572P~-57 A (inMeV) (1) 

where p is the nuclear density, p — 0.17 fm~ 3 , and f w ~ 93 MeV is the 
pion decay constant. This attraction is doubled, roughly, within chirally 
based coupled-channel KN-ttH-ttA calculations which provide also for a 
strong absorptivity Q. Shallower potentials, Re Vk-(po) ~ —(40-60) MeV 
at threshold, are obtained by requiring that the in-medium K~N t(p) ma- 
trix is derived self-consistently with the potential V K - = t(p)p it generates 
[1, 0|. In contrast, comprehensive global fits to i^~-atom strong-interaction 
shifts and widths yield very deep density dependent K~ nuclear potentials 
at threshold, in the range Re V K -(p ) ~ -(150-200) MeV Q- In this Letter 
we discuss recent FINUDA measurements that might bear on this issue by 
providing constraints on how deep Re Vk- is at threshold. 

In the preceding Letter [8;] , the FINUDA Collaboration at DA$NE, Fras- 
cati, reported on A-hypernuclear excitation spectra taken in the K~ top + A Z — > 
7r~ + j^L reaction on several p-shell nuclear targets. Formation rates were 
given per stopped K~ for bound states and for low lying continuum states. 
In the bound state formation rates agree nicely with a previous KEK 
measurement js)]. The recent FINUDA data allow for the first time to con- 
sider the A dependence of the formation rates in detail within the nuclear p 
shell where nuclear structure effects may be reliably separated out. It is our 
purpose in this companion Letter to apply one's knowledge of the nuclear 
structure aspect of the problem in order to extract the dynamical contents 
of the measured formation rates, particularly that part which concerns the 
K~ nuclear dynamics at threshold. In doing so we transform the partial 
formation rates reported for well defined and spectroscopically sound final A 
hypernuclear states into Isa hypernuclear formation rates that allow direct 
comparison with DWIA calculations. 

The expression for the formation rate of hypernuclear final state / in 
capture at rest on target g.s. i, apart from kinematical factors, is a product of 
two dynamical factors [aliol fllL Il2j: (i) the branching ratio for K~n — > n~A 
in K~ absorption at rest in the nuclear medium, here denoted BR; and (ii) 
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the absolute value squared of a DWIA amplitude given by 




(2) 



divided for a proper normalization by the integral p of the K atomic density 
overlap with the nuclear density p(r) 



Here p/, stands for the nuclear to hypernuclear transition form factor, xk s is 
an outgoing pion distorted wave generated by a pion optical potential fitted to 
scattering data, and ^f n LM is a K~ atomic wavefunction obtained by solving 
the Klein-Gordon equation with a K~ nuclear strong interaction potential 
Vk- added to the appropriate Coulomb potential. The integration on the 
r.h.s. of Eq. is confined by the bound-state form factor pfi to within 
the nucleus, where ^ n LM is primarily determined by the strong-interaction 
Vk-, although *$? n LM is an atomic wavefunction that peaks far outside the 
nucleus. The sensitivity of the DWIA amplitude Eq. (J2J) to Vk- arises from 
the interference of ^ h lm with the pion oscillatory distorted wave Xq/ • 111 
particular, once Vk- is sufficiently deep to provide a strong-interaction bound 
state for a given L, the atomic ^ u lm also becomes oscillatory within the 
nucleus which mag nifies the effects of interference, as verified in past DWIA 

calculations 0,Ej. 




In this Letter we point out another strong sensitivity to the initial-state 
K~ nuclear dynamics arising from the energy and density dependence of 
the K~n — > 7r~A BR. We show how to incorporate this energy and density 
dependence into the calculation of a properly averaged value BR which de- 
pends on the K~ atomic orbit through L and on the mass number A of the 
target. The resulting calculated Isa formation rates are then compared to 
those derived from the FINUDA data and conclusions are made on the deep 
vs. shallow K~ nuclear potential issue. 

2. Derivation of ls A capture rates from FINUDA data 

The FINUDA spectra show distinct peaks for several Is a and lp\ states 
in the nuclear p shell. In general, the derivation of the lp\ formation rate is 
ambiguous given that the lp\ formation strength is often obscured by a ris- 
ing A continuum. In A Be and in X A C it is also mixed with a substantial part 
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of the Isa formation strength owing particularly to high lying T = 1 parent 
states in the corresponding core nuclei. For this reason, we here deal only 
with the Isa formation strength, deriving it in each p-shell A hypernucleus 



from unambiguously identified low lying Is a states. According to Ref. [13 
the corresponding hypernuclear formation rates are given by a ls\ formation 
rate R(ls\), which is independent of the particular hypernuclear excitation 
considered, times a structure fraction derived from neutron pick-up spec- 
troscopic factors in the target nucleus. This theoretical framework is also 
applicable to forward cross sections of in-flight reactions such as (tt + ,K + ) 
and (e,e'K + ). In Tabled] we present Isa formation rates derived from the 
FINUDA K~ capture at rest hypernuclear spectra js), HH for a procedure 
denoted (a). In each spectrum we focus on the strongest low- lying particle- 
stable hypernuclear excitation which is also well described in terms of a A 
hyperon weakly coupled to a nuclear core parent state. These core parent 
states are listed in the table. The measured formation rates for the cor- 
responding hypernuclear excitations from Refs. jil, [l4| are then divided by 
the structure fractions listed in the table to obtain values of -R(Isa)- For 
comparison, we display in the last column the Isa component of forward- 
angle integrated (tt + , K + ) cross sections, also derived using the peaks listed 
in the second and third columns. These (tt + ,K + ) strengths show little A 
dependence, in contrast to the K~ capture at rest Isa formation rates that 
decrease by a factor 3.5 in going from 7 Li to 16 0. 



Table 1: Isa formation rates R(ls\) per stopped K , derived from the strongest hypernu- 
clear bound state peak for each of the listed targets [procedure (a)]. Data are taken from 



the preceding Letter [8(, and for * A C from [14J. The errors are statistical and systematic 



in this order. The Isa structure fractions are from [l5| and, if unlisted there, from [la ]. 
Listed in the last column, for comparison, are Isa forward-angle integrated (tt + , K + ) cross 
sections, also derived by using procedure (a) from KEK-E336 measurements [16J. 



target 


peak 


E* 

core 


ls A 


R(ls A ) x 10 3 




A Z 


core 


(MeV) 


frac. 


per stopped K~ 


(ti+,K + ) 


7 Li 


3+ 


2.19 


0.311 


1.48 ±0.16 ±0.19 


1.56 ±0.10 


9 Be 


2+ 


2.94 


0.242 


0.87 ±0.08 ±0.12 


1.40 ±0.05 


12 C 


(3/2)~ 


g.s. 


0.810 


1.25 ±0.14 ±0.12 


1.78 ±0.04 


13 C 


2+ 


4.44 


0.224 


0.85 ±0.09 ±0.13 


1.87 ±0.09 


16Q 


(3/2)- 


6.18 


0.618 


0.42 ±0.06 ±0.06 


1.47 ±0.05 
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In the second procedure, denoted (b) and presented in Table [2J we con- 
sider all the particle-stable ls\ states corresponding to observed peaks for 
which the shell model offers reliable identification. For three of the five tar- 
gets listed, this procedure saturates or is close to saturating the ls\ formation 
strength. However, in both A Be and l \C the Is a particle stable hypernuclear 
states represent less than half of the full ls\ strength. In the last column of 
Table [2] we assembled Isa forward-angle integrated (tt + ,K + ) cross sections, 
derived this time by applying procedure (b). Similarly to Table [IJ the weak 
A dependence of these Is a (tt + , K + ) cross sections is in stark contrast to the 
fast decrease of the ls\ formation rates, again by a factor 3.5, going from 7 Li 
to 16 in K~ capture at rest. The strong A dependence of the (K~ top , n~) 
rates with respect to the weak A dependence of the (tt + ,K + ) cross sections 
reflects the sizable difference between the strongly attractive K~ nuclear 
interaction at threshold and the weakly repulsive K + nuclear interaction. 

Table 2: Same as in Table [T] except for using several (rather than one) well defined Is a 
bound states for each of the listed targets [procedure (b)]. 



target 


peaks 


ISA 


R(ls A ) x 10 3 


°n A (l-ib) 


A Z 




frac. 


per stopped K~ 


(ti+,K + ) 


r Li 


1,2,3 


0.833 


1.25 ±0.14 ±0.17 


1.29 ±0.12 


9 Be 


1,2 


0.435 


0.85 ±0.09 ±0.11 


1.20 ±0.05 


12 C 


1,2,3 


0.995 


1.67 ±0.23 ±0.23 


1.92 ±0.07 


13 c 


1,2 


0.347 


0.84 ±0.12 ±0.12 


1.93 ±0.12 


16Q 


1,2 


1.000 


0.36 ±0.06 ±0.05 


1.32 ±0.05 



It is encouraging to see that both sets of R(1sa) values in Tables [H 
and [2] are consistent within statistical uncertainties with each other, except 
marginally for 12 C which dates back to a separate FINUDA run MJ. Pro- 
cedure (a) yields a value for 12 C that compares well with .R(1sa)[ 12 C] = 
(1.11 ±0.14) x 10~ 3 per -fQ™, the latter value corresponding to E x < 7 MeV 
in the KEK X A C spectrumjj Therefore, in the present study we adopt the 
.R(Isa) values listed in Tabled] 



1 We thank Dr. Tamura for providing details from his Ph.D. thesis [17j on the KEK 
experiment [9|. 
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3. Energy and density dependent K n->7r A branching ratios 




-100 -80 -60 -40 -20 -100 -80 -60 -40 -20 

E-E lh (MbV) E-E, h (MeV) 

Figure 1: Subthreshold energy dependence of the K~n — > ir~ A branching ratio BR in 
the CS30 version of the chirally motivated model Ref. 18]. The l.h.s. curves for 50, 100% 
nuclear matter density demonstrate Pauli blocking effects whereas the r.h.s. curves account 
additionally for self energy effects. 



Figure [T] shows the subthreshold energy dependence of the free-space 
K~n — > 7r~A BR generated by the CS30 version of the chirally motivated 
coupled channel model of Ref. |l8j]H This 1 = 1 BR is about 10% at thresh- 
old, decreasing to roughly half of its value as the 1 = A (1405) subthreshold 
resonance is traversed, and then increases to approximately 40% on approach- 
ing the 7rS threshold about 100 MeV below the KN threshold. The figure 
also shows the in-medium BR below threshold at densities 50% and 100% of 
nuclear matter density po, in two versions of medium modifications. In the 
l.h.s. plots the only medium effect is Pauli blocking, which acts in interme- 
diate KN states in the coupled channel equations. This is known to have 
the effect of p ushing the dynamically generated A(1405) to energies above 



threshold [19l . |20j, thus weakening the 1 = interaction and consequently 
increasing the 1=1 BR. The energy dependence in the subthreshold re- 
gion is seen to be monotonic. The r.h.s. plots include in addition to Pauli 
blocking also meson and baryon self-energy (SE) terms in intermediate state 
propagators. This pushes back the A(1405) [a, Hlj and in the chirally based 



2 The parameters of CS30 are constrained by a^N = 30 MeV. 
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model used here [18| results in stronger energy and density dependencies. 



The same chiral model was used in Ref. [12[ to generate a K~n — > n~A BR 
which, however, was fixed at its threshold value, thus neglecting any possible 
energy dependence. Since the in-medium BRs plotted in Fig. [I] exhibit a 
sizable energy and density dependence, it is essential to consider the implied 
effects in the evaluation of the Isa formation rates. 

The K~n — > n~A BR depends on the initial K~n invariant energy -y/i, 
with s = (E K + E N ) 2 — (pk+Pn) 2 in obvious notation. In the two-body cm. 
system p K + p N = 0, but in the i^~-nucleus cm. system (approximately 
nuclear lab system) p K + p N ^ and averaging over angles yields (j>k + 
Pn) 2 — > (Pa'+Ptv)- For bound hadrons, with Ek = mx—BK, En = ttin—Bn, 
we expand near threshold, E t ^ = nix + m^, neglecting quadratic terms in 
the binding energies B K , B N : 

r* ~ f r r mN mK v2r (A\ 

Vs « ±/ t h - B N - B K ■ ■ . (4J 

m N + m,K 2rriN m N + m K Lvtik 

For K~ capture at rest, we further neglect the atomic B K with respect to B N 
and replace the K~ kinetic energy p\j [2mK) in the local density approxima- 
tion by —Re Vk- (p) which dominates over the K~ Coulomb potential within 
the range of densities of interest. The neutron kinetic energy p N /(2rriN) is 
approximated in the Fermi gas model by 23 (p/po) 2 ^ 3 MeV. Altogether the 
energy argument of the K~n — > tt~A BR assumes the forrrJl 

w E th - B N - 15.1 (p/p ) 2/3 + 0.345 Re V K -(p) (in MeV) (5) 

which unambiguously prescribes the subthreshold two-body energy as a func- 
tion of nuclear density at which BR(y / i, p) of Fig. CD is to be evaluated]^ Note 
that Eq. (jSJ) leads to implicit density dependence of BR( v /s, p) through the 
invariant energy variable s/s, in addition to the explicit p dependence. The 
input BRs for our ls\ hypernuclear formation rates calculation were obtained 
by averaging this chiral-model BR(y / i, p), for a given V^-(p), over the K~ 
nuclear density overlap p(r) | \I/ n LM( r ) 1 2 of Eq. (j3J): 

BR = I ! BR(v^,p) p(r) \ * nLM (r) | 2 d 3 r. (6) 



3 Applications of this form to kaonic atoms will be discussed elsewhere [221 ] . 

4 Related ideas on the relevance of extrapolating to subthreshold energies in K ~ capture 



at rest have been repeatedly made by Wycech, see Ref. [2c 
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For B N we used target neutron separation energies. The nuclear densities 
used were obtained from modified harmonic oscillator nuclear charge densities 
by unfolding the finite size of the proton. The structure of Eqs. (jSJ), (EJ), 
together with the plots of Fig. (TJ imply that deep K ~ nuclear potentials lead 



to significantly higher values of BR than the threshold value used in Ref. [12 
which indeed is borne out by the present calculations. 

4. Confronting data with calculations 

The Is a formation rates for a shallow K~ nuclear potential V^? of depth 
— ReV^?(p = po) ~ 50 MeV and for a deep K~ nuclear potential V^? of 
depth — ReV^?(p = po) ~ 190 MeV have been recalculated with refined 
K~ atomic wavefunctions and n~ distorted wavesJl A major change here 



with respect to Refs. [12j, [24j is the use of energy and density dependent BRs 
as outlined in Sect. |3j The resulting BRs for the deep K~ potential V^P 
display considerable A dependence, with values higher than the threshold 



value used in Ref. [12|, particularly from C on. In contrast, the BRs for 
the shallow potential V^? show little A dependence, with values lower than 
the threshold value. Consequently, the difference between the DD and SH 



rates is no longer as large as calculated for a fixed BR threshold value [12 
For example, the calculated rates are (15-18)% of the experimentally derived 
rate for 7 Li under procedure (a) using the best-fit pion optical potential n e , 
Eqs. (19), (20) of Ref. [11, and (23-26)% of it using the pion optical potential 
7Tfe employed in Ref. 12| . 

Here we focus on the A dependence of the ls\ formation rates. For given 
K~ and tt~ potentials the calculated rates are scaled up by a normalization 
factor to achieve agreement for 7 Li with the Isa rate derived from the data 
under procedure (a) in Table [H This is shown in Fig. [2] where the uncer- 
tainties of the experimentally derived ls\ rates consist only of statistical 
errors that vary from one target to another. The systematic errors, on the 
other hand, are the same for all targets and drop out when considering A 
dependence within the present set of FINUDA data. The normalized calcu- 
lated Isa rates shown in the figure are based on BRs calculated according to 
Eq. © from the BRs plotted on the l.h.s. of Fig. □ (CS30, no-SE). Results 



5 The complex K nuclear potentials V^l and VBP were denoted K x and -STdd, re- 
spectively, in Ref. 113] where a complete listing of their parametrization is available. 
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Figure 2: Comparison between Isa formation rates derived from the FINUDA K~ capture 
at rest data [H, [3] and DWIA calculations normalized to the Isa formation rate of A Li 
listed in Table [1] for shallow (SH, solid) and deep (DD, dashed) K~ nuclear potentials. 
The calculated Isa formation rates use K~n — > ir~ A in-medium BRs without self energies, 
see Sect. [31 and pion optical potential 7r e from Ref. [f|. The error bars consist of statistical 
uncertainties only. 



are shown for the pion opticalpotential ir e which was fitted to 7r~- 12 C angu- 
lar distributions at 162 MeV [6j, and for the two K~ nuclear potentials V^? 
and V^?P. We note that the decrease of the experimentally derived Is a rates 
from 7 Li to 9 Be, followed by increase for 12 C and subsequently decreasing 
through 13 C down to 16 0, is well reproduced by both calculations shown in 
Fig. [2J However, the deep V^?P calculated rates reproduce better the A de- 
pendence of the experimentally derived rates than the shallow V£r potential 
does. In reaching this conclusion on Vj^P, the increase of the K~n — > tt~A 
BR values between 7 Li and 16 is essential, by moderating the fall off of the 
rates calculated using A independent BRs. Similar conclusions hold for the 
A dependence of rates calculated using BRs that are derived according to 
Eq. © from the BRs plotted on the r.h.s. of Fig. Q] (CS30, with SE). On 
the other hand, if the pion optical potentials 7r& or n c (applied in Ref. [l2j]) 
are used in these calculations, neither VPP nor V^P do as good a job as the 
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combination V^P and the best-fit ir e does, and no firm conclusion can be 
drawn. 



5. Conclusion 

In conclusion, we have derived Isa hypernuclear formation rates from 
peak formation rates associated with the A Z(K~ t , vr~)^Z spectra presented 
recently by the FINUDA Collaboration on several nuclear targets in the 
p shell jsl, [l4| • We then compared the A dependence of these derived rates 
with that provided by calculations for the two extreme Vk- scenarios dis- 
cussed at present, a shallow potential 0, @] and a density dependent deep 
potential [7J. The calculations use K~n — > n~A in-medium BRs generated 
by applying a recent chirally motivated coupled channel model [18|. These 
BRs exhibit a strong subthreshold energy and density dependence, as shown 
in Fig. [H and therefore result in A dependent input values BR that depend 
sensitively on the initial-state K~ nuclear potential Vk-- The calculations 
also demonstrate additional strong sensitivity to Vk- through the atomic 
wavefunctions it generates which enter the DWIA amplitude Eq. (J2J), as dis- 
cussed extensively in previous calculations of K~ t hypernuclear formation 
rates [6, 12|. The comparison between the calculated A dependence and that 



derived from the FINUDA data slightly favors a deep K~ nuclear potential 
over a shallow one. This conclusion outdates the one reached in an earlier 
version in which the energy and density dependence of the BRs, resulting 
here in a new source of sensitivity to Vk--, was disregarded [241] . In future 
work, it would be interesting to use other versions of K~N chirally moti- 
vated models and to extend the range of nuclear targets used in stopped K~ 
reactions to medium and heavy weight nuclei in order to confirm the present 
conclusion, and to look for more subtle effects of density dependence. 
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